The Gribov horizon and spontaneous BRST symmetry breaking 



David Dudal* 

Ghent University, Department of Physics and Astronomy, Krijgslaan 281-S9, 9000 Gent, Belgium 

Silvio Paolo Sorella^ 

UERJ - Universidade do Estado do Rio de Janeiro, Instituto de Fisica, Departamento de Fi'sica Teorica, 
Rua Sao Francisco Xavier 524, 20550-013 Maracana, Rio de Janeiro, RJ, Brasil 

An equivalent formulation of the Gribov-Zwanziger theory accounting for the gauge fixing ambiguity in the 
Landau gauge is presented. The resulting action is constrained by a Slavnov-Taylor identity stemming from a 
nilpotent exact BRST invariance which is spontaneously broken due to the presence of the Gribov horizon. This 
spontaneous symmetry breaking can be described in a purely algebraic way through the introduction of a pair 
of auxiliary fields which give rise to a set of linearly broken Ward identities. The Goldstone sector turns out 
to be decoupled. The underlying exact nilpotent BRST invariance allows to employ BRST cohomology tools 
within the Gribov horizon to identify renormalizable extensions of gauge invariant operators. Using a simple 
toy model and appropriate Dirac bracket quantization, we discuss the time-evolution invariance of the operator 
cohomology. We further comment on the unitarity issue in a confining theory, and stress that BRST cohomology 
■ alone is not sufficient to ensure unitarity, a fact, although well known, frequently ignored. 

: 

(N . 

->-> ' I. INTRODUCTION 

<D ■ 

The issue of the BRST symmetry in presence of the Gribov horizon is of pivotal relevance in order to achieve a better 
l— ~~ 1 1 understanding of the infrared nonperturbative region of a confining and ultraviolet asymptotically free Yang-Mills theory. At the 
moment, it seems safe to state that reconciling BRST symmetry with the Gribov horizon is a task not yet fully accomplished. As 
one learns from the original works of Gribov and Zwanziger [1-3], the restriction of the domain of integration in the functional 
Euclidean integral to the first Gribov horizon enables one to take into account the existence of the Gribov gauge fixing ambiguities 
in a local and renormalizable framework, as summarized by the Gribov-Zwanziger action in the Landau gauge [2, 3] and its more 
Q\ , recent Refined version [4-6] . 
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Two important physical consequences which arise from the restriction to the Gribov horizon are: the appearance of a nonper- 
turbative dynamical mass parameter y 2 ~ A^ CD , called the Gribov mass, and a geometrical picture for the confinement of gluons. 
Due to the presence of the Gribov horizon, the two-point gluon correlation function has no particle interpretation, displaying a 
violation of reflection positivity, a feature confirmed by numerical simulations [7]. It is worth mentioning that the expression for 
the tree level gluon propagator obtained from the Refined Gribov-Zwanziger (RGZ) theory [4-6] is in very good agreement with 
the most recent lattice data on very large lattices [8-1 1]. This propagator exhibits complex poles in momentum space 1 , a feature 
which has a natural accommodation in terms of gluon confinement and absence of corresponding asymptotic states. Despite the 
appearance of complex poles, the aforementioned confining gluon propagator has been employed to investigate the correlation 
functions of gauge invariant composite glueball operators in order to get estimates of the glueball masses. It turns out that a 
Kallen-Lehmann spectral representation with positive spectral density can be extracted from these correlation functions [12]. 
The resulting mass estimates for the lowest glueball states, ++ ,2 ++ ,0 h , are in agreement with the available numerical data 
on the spectrum of the glueballs [13]. It is worth mentioning that such type of gluon propagator has also been used in previous 
studies in hadron physics [14]. 

Although the above mentioned results can be interpreted as encouraging evidence of the fact that the Refined Gribov- 
Zwanziger theory can be effectively employed to investigate the physical spectrum of a confining Yang-Mills theory, we are 
far from having at our disposal a systematic way of constructing a set of physical composite operators whose correlation func- 
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1 To avoid confusion, let us mention here that upon including perturbative corrections, the properties of the RGZ gluon propagator in the complex momentum 
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tions can be directly linked with the observable spectrum of the theory. In particular, the characterization of the spectral analytical 
properties and of the positivity of the corresponding correlation functions appears to be a highly cumbersome problem, taking 
into account that explicit calculations have to be done by employing a confining gluon propagator exhibiting complex poles. 
Let us point out here that also other continuum approaches can adequately reproduce the lattice results for the gluon propagator 
[15-19], however the resulting correlation functions will also exhibit a complicated structure, so it is in any case unclear how 
these nontrivial propagators and eventually also vertex functions will conspire to form e.g. a truly physical glueball correlation 
function. 

We expect that a better understanding of the role of the BRST symmetry in confining gauge theories could shed some light on 
these issues. In the present paper we aim at presenting new results on the BRST symmetry in presence of the Gribov horizon 2 . 
For the benefit of the reader, let us proceed by providing a short summary of previous attempts on this topic. We may start by 
remarking that already in the original paper by Zwanziger on the renormalizability of the Gribov-Zwanziger theory [2, 3], it was 
pointed out that the restriction to the Gribov horizon induces a soft breaking of the BRST symmetry which does not jeopardize the 
renormalizability of the Gribov-Zwanziger action. Later on, in [24], an attempt was made to argue that the Gribov-Zwanziger 
action displays an exact BRST invariance which is, however, spontaneously broken. Moreover, the spontaneous symmetry 
breaking was encoded into the BRST transformations through the introduction of terms containing an explicit dependence from 
the space-time coordinates x M , see [25] for more details. As pointed out in [5], this kind of .^-dependent transformations 
is plagued by several drawbacks as, for example: the explicit loss of translational invariance and the appearance of nontrivial 
surface terms requiring appropriate boundary conditions for the fields at infinity. In [26], the soft breaking induced by the Gribov 
horizon was discussed in terms of the BRST cohomology. It was shown that the soft breaking enables to give a cohomological 
proof of the fact that the Gribov mass y 2 is a physical relevant parameter. Also, the idea of the soft BRST breaking was employed 
to construct a renormalizable quark model [27] accounting for the lattice data on the quark propagator [28, 29]. Further, in [30] 
it was pointed out that the soft breaking induced by the Gribov horizon could be converted into an exact invariance which, 
however, gives rise to a non-local and non-nilpotent BRST operator. This non-local formulation of the BRST invariance was 
subsequently casted in a local framework in [31]. Though, the resulting BRST operator still lacks nilpotency. We mention that 
in [32] a connection between the soft BRST breaking appearing in the Gribov-Zwanziger theory and the Kugo-Ojima criterion 
[33] was established. More recently, the Gribov-Zwanziger theory was reformulated in such a way that the soft BRST breaking 
was transformed into a linear breaking [34, 35], a feature which allowed to clarify the nonrenormalizability properties of the 
Gribov parameter y 2 . 

In this work we shall present a new equivalent formulation of the Gribov-Zwanziger action which exhibits an exact BRST 
invariance, with a nilpotent BRST operator. Moreover, the BRST symmetry is spontaneously broken, the spontaneous break- 
ing parameter being nothing but the Gribov mass y 2 . However, unlike the previous formulation of [24], no terms with explicit 
x^-dependence are introduced in the BRST transformations of the fields. This is an important feature of this new formulation, al- 
lowing for a fully local, covariant and renormalizable framework. Another important outcome is the decoupling of the Goldstone 
mode associated to the spontaneous symmetry breaking. Moreover, the nilpotency of the BRST symmetry enables us to make 
use of cohomology tools to identify a set of local and renormalizable operators in presence of the Gribov horizon. The paper is 
organized as follows. Section II is devoted to the presentation of the new equivalent formulation of the Gribov-Zwanziger action 
and of the spontaneous symmetry breaking of the BRST operator. In Section III we collect the set of Ward identities fulfilled 
by the reformulated action and we analyze the decoupling of the Goldstone mode. In Section IV we discuss, with the help of a 
toy model, the invariance under time evolution of the set of composite operators corresponding to the cohomology classes of the 
BRST operator. Section V collects our conclusion. 



2 In the very recent paper [20], it is mentioned that a Landau gauge has been constructed which involves averaging over all Gribov copies. The latter method- 
ology would be an alternative to handling the Gribov copy problem compared to the restriction within the Gribov horizon. Again according to [20], it has 
been shown in some of the quoted literature that such averaging harbors an intact BRST symmetry. We have however been unable to recover from that quoted 
literature either how to average over all Gribov copies in the continuum, or why this would induce a well-defined BRST symmetry. In fact, a recent continuum 
attempt towards sampling over all Gribov copies in [21] seems to give rise to a Lagrangian without BRST invariance, given the emerging connection with the 
Curci-Ferrari model, which is known to not exhibit exact BRST invariance. A similar observation was made in [22]: "Only if the numerical lattice procedure 
converges towards a correct sampling of the fundamental modular region in the infinite volume and continuum limits, can the BRST breaking effects be 
expected to go away.". Clearly, the disappearance of BRST breaking effects in the continuum is to date thus merely a wishful expectation, being an extremely 
hard problem to actually solve. No continuum approach to restricting the functional integration to this fundamental modular region [23] is known. 
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II. REFORMULATION OF THE GRIBOV-ZWANZIGER ACTION AND SPONTANEOUS BRST SYMMETRY BREAKING 

A. An equivalent reformulation of the Gribov-Zwanziger action 

Our starting point is the Gribov-Zwanziger action which enables us to restrict the Euclidean functional integral to the first 
Gribov horizon, namely 

Sgz = \\ A W + / d rf JC (fd^ + ^Df c b ) + J d^x (^'3^7(0^ 

+ J d'x^gf^A^+^l-diN 2 -!)^) . (1) 



= 0, 



The Gribov parameter y 2 is not a free mass, but it is determined in a self-consistent way through the gap equation -J^f 
where Zs vac stands for the vacuum energy of the theory [1-3]. 

The foregoing model has a softly broken BRST invariance, generated by the nilpotent transformations 

sAl = -Dfc J ' = -(d fl 8 ab +gf acb A c ^c b , sc a = %-f ach c h c c , sc a = b a , sb a = , 

s< = qjf , - , stf = < , ,< = . (2) 

Indeed, we have s 2 — while 

sSgz = -T 2 / & d x(gf ahc Dfc\t; + ^+gf abc A a ^ b ;) . (3) 
We propose to replace expression (1) by the following action 

+ J d d x (y a ;dvD a v m <?; K - m» c d v D™(»z c - 8 (d v o%) f abm (D vC ) b (p- c ) 

/, / -~~ a b „ ^ , ab ^ , ab , , , -~~ab ... , . Cl ab \ 

& d x (- q^ 2 g; b + t^t* + r^Dftof g^Df^ + gf M f^ c y^ 

+ J d d x (g% s.vS-Y) + (C - 5,v8" Y) - sit*) (4) 



by introducing two new BRST quartets of fields 3 consisting of Q^ b , 7^ (commuting), jfjffi, (anticommuting) and their 
hat-counterparts: 



7r ah 


-p^ab 
- ^uv> 




= o, 


C: ab 


^ab 

= G/jv' 


^ab 


= o, 






_^ crab 


v arab 


= o, 


crab 


_ rrab 

J fjV 1 


s 7rab 


= 


(5) 



and the singlets Mffi , with thus s!Hffi = sJ^ b = 0. As the action (4) can be easily rewritten as 
Sgz = \j A d xF^ v + j d'xs^d^ + ^D^) 

+ j d d xs {-r^Qt -f&3# -K d ?<) 

+ j d d x - s.vS^Y) + (C - M fl V) - ^vC) . (6) 



3 We dressed these new fields with the designated color and Lorentz-stracture as this formulation allows for several Ward identities, which will be employed to 
prove the all orders renormalizability of the action S' az in a follow-up paper. 
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the BRST invariance is almost immediate, sS' GZ = 0, whereby we have preserved the nilpotency of the BRST operator, s 2 = 0. 

-at ^- ab 
J/JV ~ Q 'ijV ~ 



The role of the fields is clearly seen, as upon using their algebraic exact equations of motion (EOM), Q = Q 



y 2 8 a/ '8 jU v, we get as an immediate consequence 

' | d d x (-f g f bc A; (<# + <p£) - d(N 2 -iyf). (?) 



Naively, the doublets (jF, (7), (jF, (7) together with the i#)-terms would be sufficient to write down a BRST-invariant action 
that reduces, upon using the (9{,y{)-EOM, to the original GZ action (1). However, this would leave us with a trivial partition 
function, due to the vanishing path integration over [&7] or [djF]. The use of the double quartet circumvents this problem, 
neither corrupting the BRST invariance nor the exact correspondence with the action (1). In fact, the integration over f , J , f 

and f will now generate a unity, due to absence of any other term in f or thereby removing all terms ~ f and ~ f . This 
can be best seen by performing 3 simultaneous shifts of integration variables: 

7$ -> ?ft - ^ {^ b +gf ak % p c^ b ) , < -+ < + [m- l f (gf aU gp> e M p c p ) , 

rfi -> ^ - ^ (rf<tf+gf M i%<?$?) , (8) 

from which the equivalence between the two formulations, (1) and (4), becomes apparent. 



B. Spontaneous symmetry breaking of the BRST symmetry 



In order to realize that the new formulation of the Gribov-Zwanziger action does exhibit a spontaneous symmetry breaking of 
the BRST symmetry, we rewrite the action S' az by making explicit use of the equations of motion of the auxiliary fields {M, 9€): 



~C = Q a l = i^- (9) 



Thus 



+ J d d x (jja 2 ^-? 2 . (io) 

Expression (10) is easily seen to be left invariant by the following nilpotent BRST transformations: 

sAl = -Dfc b = -(d fl 5 ab + gf acb A'j J )c b , sc a = ^f ach c h c c , sc" = b a , sb a = , 



sWf = cpf , s% b - , ,cpf = < , srf = , 

-_ab 



= , s?"^ = fS?*8v , sQ ab = 7$ , s7$ = 0, sQ ab = <F« b , = (11) 

with 

sS' GZ = 0, s 2 = 0. (12) 
Moreover, from equations (1 1), it follows that the BRST operator suffers from spontaneous symmetry breaking, namely 

= Y^V , (sf%) = Y>S a %v ■ (13) 

A few remarks are in order: 
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• unlike the formulation of [24], terms with explicit dependence from the space-time coordinates x^ are introduced neither in 
the BRST transformations (11) nor in eqs.(13). This relevant feature enables us to study the spontaneous breaking within 
the well-established framework of genuine local quantum field theory. Indeed, one notices that both the action S' GZ , eq.(4), 
and the BRST transformations, eqs.(l 1), are local. 

• as it is apparent from eqs.(13), the spontaneous breaking of the BRST symmetry is entirely driven by the Gribov mass 
parameter y 2 . This means that the breaking is due to the restriction to the Gribov horizon. As such, it is a truly non- 
perturbative effect with a clear geometrical meaning. 

• it is worth to point out that the introduction of the auxiliary fields (?{, H) enables us to treat the spontaneous symmetry 
breaking in a purely algebraic way, as it follows from the equations of motion of these fields: 

S.vS-Y , 5k = C - M"Y • (14) 

In fact, according to the algebraic renormalization procedure [36], equations (14) will acquire the meaning of linearly bro- 
ken Ward identities at quantum level, providing a very powerful and elegant treatment of the properties of the spontaneous 
symmetry breaking at the quantum level, see next Section, 

• although not worked out in explicit way here, we underline that the spontaneous symmetry breaking mechanism presented 
in this work generalizes without any difficulty to the Refined Gribov-Zwanziger action [4-6]. Also a soft breaking in the 
fermion sector [26, 27] fits within this new scheme. 



III. RESTRICTIONS ON THE EFFECTIVE ACTION 



A. Set of Ward identities 



In a follow-up paper [37], we will discuss the renormalizability of the new reformulated GZ action (4). For that purpose, let 
us already list here a set of Ward identities that will be helpful to construct the most general admissible counterterm. According 
to the algebraic renormalization setup [36], we have to couple the nonlinear BRST variation of and c a to BRST invariant 
sources K£ and L". In addition, following [2, 3], we can introduce a composite index notation, ;' = (v,b), as follows 

%t = Tl? = Ttf, (15) 

where X^ b ,,"(\" h stand for generic fields of the action (4). Also, for the identification of two more powerful Ward identities, we 
have found it useful to introduce a further BRST doublet of external sources 

sXf = pf, spf = (16) 

coupled to the composite operator f^D^cr and its BRST variation Q fji D bk c k . 
We thus look at Ward identities for the complete classical BRST invariant action 

L = ly d d xF^F« v + J d d x {b a d^ + ^Dfc b ) 

+ j d d x (V? 3 v D™(p™ - C5f dyD^a/p - g (d v &f)f abm D b v k c k ^p) 

+ J d d x (-f > 2 + f 1^9" + %Df(o k %ptf$ + gf ak( %D ( ; c ^ k ) 

+ J d d x ( 0» - 8;, y 2 ) + % (|; - S^.y 2 ) - ^ ( .f ;.) 

+ J d d x (-K«Dffc b + Uff ahc c h c c + pf^\p bk c k - XffliD bk c k ^ . (17) 
It turns out that this action E fulfills the following Ward identities: 
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1. An exact Slavnov -Taylor identity 

r d I 8E 8E 8E 8E 8E _^ 8E a 8E —a 8E 8E ^> 8E ~ 8E ab 8 

7 ^8A« + 8Z»8^ + 8^ + (p( '8sf + G>i 5^f + ^"SjJ + " i 8^ + "'$f\ + "slft"^' 



= 0. 
(18) 



2. 77ze linearly broken Ward identity for the Gribov mass parameter J 2 

8E 



W = ~fd 4 x(?Q + ^). (19) 



3. The linearly broken equations-of-motion identities 



g-W g+^=0, (20) 
8E 2 — a 8E 



5?;, " s ^ = " 8 ? "" <21) 

£--S»T„ Jr=^. < 22 > 

°-V 8^,- Sj7> 

5X Sj-tffc, -5-=f£-tfft. (23) 



8^" 



4. Other linearly broken equations-of-motion type identities 

8W +3v 8^. Dv 8^. 



^+3 v -^-Df^ = -d^ + ^K-d^-gf^Ai^f, (24) 



8E 8E w / 8E 



80? - ^ 

5. 77ze linearly broken integrated <$-equation-of-motion 



6. 77ie exacf integrated Ward identities, generated by: 



8 


^ 8 






8 


-pa 8 


8 % 


w 8 J fl 



8 ^."8 



- (/„— I * = o. ( - E = 0, 

rfl 8 8 



7. 77ze linearly broken integrated global Q { j Ward identity 

8 8 „ 8 _„ 8 — fl 8 
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8. An integrated Ward identity involving the Faddeev-Popov ghosts and the Gribov-Zwanziger fields 




All these identities are compatible with the Quantum Action Principle [36]. As such, they will severely constrain the pos- 
sible counterterm which may arise at the quantum level. Although the output of the algebraic renormalization analysis will 
be presented in a more complete and detailed work, we expect that this set of Ward identities will ensure the multiplicative 
renormalizability of the reformulated GZ action (4), as supported by the physical equivalence between the new and the original 
formulations of the GZ theory. In particular, this will imply that the quantum effective action, T = L+0(h), will fulfill the same 
set of Ward identities obeyed by £. 



B. Decoupling of BRST Goldstone mode 



The exact quantum EOMs, = = (f - - fh a -, = -4E- = Q. - f&> will lead to = Q . = fd a -, and thus to the 

spontaneous breakdown of the continuous BRST symmetry, henceforth a Goldstone degree of freedom must emerge. The latter 
can be created out from the vacuum at no energy cost using the BRST current. Writing down here only the relevant part of this 
current 

rf RST = ...+ TM a > + > (30) 

it is apparent that in the condensed vacuum, a linear combination of + enables us to identify the BRST Goldstone degree 
of freedom. However, the quantum version of the Ward identities (22) immediately lead to 

— = ^ = -< i7 a 2 8(*-y) , (31) 



telling us that the Goldstone mode is exactly a free particle. As such, it can be trivially decoupled from any consistent definition 
of a possible physical subspace. 

Let us also point out that, relying on the observations that the spontaneous symmetry braking of the BRST symmetry can 

be encoded into the (0-C, i?/)-EOMs, which have the meaning of Ward identities for algebraic degrees of freedom, and that the 
ensuing Goldstone degree of freedom itself is also completely characterized by a Ward identity, we can speak about an algebraic 
breaking of the BRST symmetry. 



IV. INVARIANCE UNDER TIME EVOLUTION OF THE BRST OPERATOR COHOMOLOGY 



A. Preliminaries 



The aim of this section is to exploit some consequences of having at our disposal an exact BRST invariance which is spon- 
taneously broken. In particular, the nilpotency of the BRST operator, eqs.(2),(5), opens the door to make use of cohomology 
tools in order to identify the BRST invariant local composite operators. From the well established results on the cohomology of 
Yang-Mills theories, see [36, 38], it follows that the cohomology classes in the space of zero ghost number are given by local 
colorless gauge invariant operators built up with the field strength F^ v and its covariant derivative D^ b . Up to this point things 
look similar to the usual unbroken case. It is worth to emphasize there is, however, a deep difference between the unbroken and 
the broken case. The BRST invariant operators are always defined modulo exact BRST terms which, in the unbroken case, are 
physically irrelevant. This is not the case in the present situation. The exact BRST pieces might now give nontrivial contributions 
to the correlation functions, due to the non-invariance of the vacuum. A detailed study of how these exact terms contribute to 
a given correlation function is beyond the scope of the present letter. Though, as a first result, we shall give here a proof of the 
invariance under time evolution of the cohomology of the BRST operator in the spontaneously broken case. Denoting by {0} 
the set of local BRST invariant operators, including the exact terms, it turns out that { 0} is closed under time evolution, namely 
the evolution in time of { 0} will not bring us outside of { 0}. 
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In order to give a proof of this statement, we shall make use of a toy model, already introduced in [12], which reproduces 
the main features of the Gribov-Zwanziger theory without leading to notational clutter with indices etc. It is specified by the 
non-local Euclidean action 

S = /d 4 *i ¥ (-3 2 -|?) ¥ , (32) 
which is easily seen to give rise to a Gribov type propagator 

M*)Y(-*)> = ■ 03) 

Despite its non-locality, the action (32) can be cast in local form, which reads 

S = yd 4 xQ\|/(-3 2 )\|/+9(-3 2 )(p + c5(a 2 )(»+eV(9-9)) (34) 

whereby we introduced 2 commuting (9,9) and 2 anti-commuting scalar fields (<0,<5), with BRST transformation 

s\|/ = 0, s(f> = co, .s(0 = 0, sa> = (f>, scp = 0. (35) 

We notice that, as in the case of the Gribov-Zwanziger action (1), the BRST symmetry is softly broken by the mass mixing 
term. Nevertheless, we can repeat the same procedure employed to construct the action (4), and reformulate expression (34) in a 
manifest BRST invariant fashion. As done before, we introduce an additional BRST quartet (jF, !F, Q, Q), and a singlet H, and 
consider 

s = ^d 4 xQ\(/(-a 2 )\(/+(p(-a 2 )(p+co(a 2 )co+e 2 \(/((p-(p) 

+g(-d 2 )g+J(d 2 )f + £(9-9)- j ¥ co+^-e 2 )) , 06) 

with 

sQ = F, sf = 0, s~J = ~Q, s~g = 0, sM = 0. (37) 

It is easily checked that sS — 0, while the i^-EOM will set Q — 8 2 , and by integrating out the (!F,!F) fields, equivalence at the 
quantum level with the original local action (34) is again established. Here, the iT-field plays the role of the Goldstone mode for 
the BRST breaking, which is linked to the exact quantum EOM, ^ = 0. The BRST current reads 

j„ = - (9„<p) oo + (3/1®)?+ @pF)G - {^G)F (38) 
Let us mention that the action (36) seems to harbor a local "gauge" symmetry, generated by 

8 £ ^(x) = 3 2 e(x) , 8 £ g(x) = e(x) (39) 
with e(x) an infinitesimal function of x. This is actually corresponding to the following local Ward identity: 

It is also there in the reformulated version of the GZ action, one just needs to combine the identities (21)-(23). In order to avoid 
to deal with first class constraints — which typically correspond to a gauge freedom [39, 40] — we shall canonically quantize the 
theory (36) without the explicit use of the ^/-field in the Hamiltonian formalism, and we shall impose Q = 2 as the primary 
constraint to start from. 



B. A few words on Euclidean time evolution 



Strictly speaking, we are not allowed to Wick-rotate a Gribov-like action (36) from Euclidean to Minkowski space. We need 
to work directly in Euclidean space. With the normal connection xo = — ix' , where the prime refers to Euclidean space, we find 
that the Euclidean Heisenberg equation, which governs the time evolution of any operator O, is given by 

grO— lO,H|, (41) 
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with solution 

O(x' ,x) = e x ° H 0(x)e- x 'o H . (42) 

Let us also have a look at the definition of the Euclidean Hamiltonian. We remind that the Euclidean action is given by 

L' = -L(... = -*'... o) (43) 

where each time component needs to be transformed from Minkowski to Euclidean signature in the designated way. In 
Minkowski space-time, one introduces the conjugate momentum of a field 4 (j) as 

8L 8L' , .3 3 

11 .0 >W iU - * = a^' * = (44) 

where we introduced the Euclidean conjugate momenta II' in a natural way from the Euclidean Lagrangian. For the Hamiltonian, 
we then find 

H = J dx (<j>n-£) = j dx {-tyU' + L') . (45) 

At the level of the quantization, the usual equal-time-commutator reads 

[(j)(x,x ),n(y,x )] = /8(x-y) . (46) 

Thus at the Euclidean level, we shall set 

[V(x,x> Q ),Yl'(y,x' )] =5(x-y). (47) 



C. A few words on the definition of the charge 



Sometimes one reads that in case of spontaneous symmetry breaking, associated to a conserved current 7^, the following 
definition of the charge , 

Q=J&\J Q , (48) 

fails. This is true, in the sense that it is not possible to define a global version of the charge, since the behavior of the massless 
modes at infinity jeopardizes the convergence of the integral. However, it should still be possible to consider the local variation 
of any quantity w.r.t. the transformation behind the conserved current 5 . As explained in Ch. 17 of [41], one defines 

(l R = Jd 4 xf R (x)a(xo)Jo(x,x Q ) (49) 

where the local action of the transformation acting on the operator is obtained from 

80= lim [Q*,0]. (50) 

R— 

In the above, one sets 

f(x)=f(\x\/R), f(y) = l if|y|<l, /(y)=0 if y > 1 (51) 

and for a(xo) one may take any smooth approximation of the 8-function, keeping / dxoa(xo) = 1. The presence of a(xo) is 
a technical issue. It amounts to a smearing in time, a feature needed for a mathematically well defined object in the case of 
relativistic QFT, in the strict Wightman-axiomatic approach to QFT We shall not be concerned with this smearing in this note 6 . 



We consider a scalar field for simplicity here. In case of e.g. a vector field, one needs to keep in mind to also transform the O-component of the field itself in 
the appropriate way. 

5 Or even more generally, to define the variation of e.g. the Hamiltonian under any transformation, symmetry-related or not. 

6 One shows that the limit R — > °° exists, and that it is independent of the choice of the smearing function a(xo) [41]. 
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We shall tacitly assume its presence and therefore simply write / d 3 x/ fl (x). The most important part of the foregoing construction 
is the presence of a cut-off at large |x| of the integrand, avoiding the potential problem of long-wavelength divergences caused 
by the massless modes. 

Thus, having introduced and the transformation law (50), we are equipped to define an operator subspace $ subject to 

G 3 «*■ [CL 0] = lim [<Z R , 0] = 0. (52) 

R — 

This subspace is then invariant under time evolution. Indeed, we easily find, using (42) and dropping from now on the prime- 
notation for Euclidean quantities, 

lim [Q/V° H 0-*° H ] = e x ° H ( lim [Q?, 0}) = 0, (53) 

as OeJ. If happens to be a nilpotent operator, the operator cohomology is invariant under time evolution. Notice that 
nowhere in this argument we used that lim^oo corresponds to a broken or unbroken symmetry generator, in the sense that 
we did not use the fact that there exists (or not) an operator SI with (limR^oofQ^,-^]) ^ 0. 

It is only at the level of taking expectation values of products of operators belonging to $ that the effects of the broken 
symmetry will emerge, as the expectation values need to be taken w.r.t the correct vacuum. For instance, as stated previously, 
in the case of the spontaneous symmetry breaking of the BRST symmetry, the choice of the nontrivial vacuum will have the 
consequence that in expressions of the form 

( [Sit (x) + sR 2 (x)] [SI, (y) + sSU (y)]> , (54) 
the s-exact pieces will give a nonvanishing contribution to the correlation function. 



D. Back to the quantization of the toy model 



Using the tools of the previous section, we first derive the conjugate momenta, 



n 7 



-doJ, 



n 



9 = 3o<p, 



I% = 3 cp, 11(0 = -d (0, Yl w = d (a, 



n. 



n 7 



(55) 



Next, we impose Q — Q 2 = as a so-called primary constraint [39, 40]. Consistency of this constraint with time evolution 7 will 
give rise to higher order "secondary" constraints. Therefore, we need the Hamiltonian, which is found to be 



+dt gd[g + djfdiW + £\|/(<P - 9) - Jv«>] 

and the Poisson bracket, which is in general given by 

{A(x,0,B(y,0} = L J d 3 z 



8(j)(z,f) 8n<^(z,f) 8n ( | ) (z,f) 80(z,f) 



(56) 
(57) 

(58) 



with t =xq identified as the Euclidean time. As we shall see, there will be no first-class constraints 8 . 

Let us now construct the secondary constraints. From Q — 9 2 = 0, we get = via { £r,H} ~ 0. Also, when Q = 8 2 , we 
have simultaneously IT^- w which, in turn, leads to d 2 Q = via |n^-,H| w 0. Here, we can stop since {8?^,H} °= d 2 H-g 
and the latter is indeed weakly zero due to a constraint already found. Likewise, {n^,H} °< —3? Q w 0, not providing any new 



information. 



7 This amounts to having a weakly vanishing (ks 0) Poisson bracket of the constraint(s) with the Hamiltonian, where weakly means modulo the previous 
constraints. 

8 First-class constraints are those with vanishing Poisson bracket with all other constraints. Second-class constraints are then the remaining ones. 
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The constraints Q — 9 2 , IT^ and 3? £7, II^ are pairwise second class, as a quick computation learns that 
I Q(x,t) — 8 2 ,n^-(y,f)| = 8(x - y) and {d? g(x,t),Hg(y,t)} = d 2 8(x — y). Any other bracket of constraints is vanishing. 

In order to quantize the theory, we first need to replace the Poisson bracket with its more general Dirac counterpart [39, 40]. 
In general, the Dirac bracket is defined as 9 



{A(x),B(y)}* = {A(x),B(y)}-£ [ d 3 z / dV{A(x),4> a (z)}fW ap 1 (z,z'){*p(z'),B(y)} . 



(59) 



Here, <I> a are the second class constraints, with the anti-symmetric invertible matrix iMap = {<t> a ,<t>|3}. In the above case, using 
the entries Q, Yl-g, djg, Ylg, this matrix and its inverse are found to be 



( S(z-z') \ 

-S(z-z') 

8 2 8(z-z') 

\ -3 2 8(z-z') / 



M- 1 





S(z-z') 


V 



-S(z-z') 









1 1 

4jt |z-z'| 






1 1 

4jt |z-z'| 



/ 



This leads to the usual brackets for all fields, with the exception of 

{^(x),r%(y)}*=0, {£(x),n ? (y)}*=0 



(60) 



(61) 



as follows after a little computation based on (59),(60). Here, we clearly see the purpose of the Dirac bracket: it becomes 
zero whenever it contains a 2nd class constraint. In order to complete the quantization of the model, we must replace the 
Dirac brackets with (anti)commutators, according to the prescription (47). Since the Dirac bracket of a constraint automatically 
vanishes, no inconsistencies will be retrieved. In fact, one has the choice to first use the constraints and then compute the bracket, 
or vice versa. Concretely, we have the following (anti-)commutators 

[<j)(x,f),n^(y,f)] ± = 8(x-y), forcj) G {\|/,(p,(p,co,co, 9" ,7} , 

[g(x,t),Ug(y,t)] = [^(x,0,%(y,0] =0. (62) 
As a last step, we need to establish that the Hamiltonian still commutes with the BRST charge, for which we rely on 

d R = J d 3 xf(x) [-n^o+n^+n^-n^] . (63) 

After some algebra, one finds 

[Q*,H] = Jd 3 xf R (x)(U Y (x,t)U g (x,t)+d l iF(x,t)d i g(x,t)] 

» 0, (64) 
upon using the constraints. The same result is obtained if the constraints are used before evaluating the bracket. 



V. CONCLUSION 



We believe the essence of the novel formulation of the Gribov-Zwanziger action presented here is that: (1) we can introduce 
the cohomology of the BRST operator which enables us to renormalize quantum extensions of the classically gauge invariant 
operators, (2) the spontaneous symmetry breaking of the BRST symmetry can be completely kept under control in an algebraic 
way, with a decoupled Goldstone mode, and (3) the operator cohomology forms a closed subspace under time-evolution, despite 
that the BRST symmetry is algebraically broken in the true vacuum. The "only" noticeable consequence of this breaking will 
be the nonvanishing of correlation functions of BRST exact quantities. One point we, willingly, did not touch upon in this letter 
is what happens in the case of spontaneously broken BRST symmetry with the gauge parameter independence of correlation 



For notational simplicity, we dropped the (equal) time t everywhere. 



12 



functions of quantities that belong to the BRST cohomology. Here, we cannot provide an answer, as the Gribov-Zwanziger 
construction is strictly limited to the Landau gauge 10 . In particular, the defining properties of the Landau gauge are heavily 
used in the derivation of the partition function. For example, a key role is played by the Hermitian character of the Faddeev- 
Popov operator, something which holds only in the case of the Landau gauge. As such, in the absence of a Gribov-Zwanziger 
construction for e.g. a general linear covariant gauge, it is meaningless trying to establish the gauge parameter independence of 
certain quantities. One could easily imagine the Gribov parameter to become gauge parameter dependent in such a setting, if 
the underlying action would be known to begin with. Nonetheless, if someday somebody would manage to construct a sensible 
Gribov-Zwanziger action with gauge parameters, the exact Slavnov-Taylor identity established here could be used to derive a 
kind of Nielsen identities [42] to control the gauge parameter dependence. 

At this point, we have reached the same conclusion as anybody can reach using the QCD action in a particular gauge which 
would have perfect nilpotent BRST invariance: we have identified a well-defined subspace of renormalizable quantum extensions 
of the classically gauge invariant operators. This is the first step in the fundamental issue of establishing unitarity of the confined 
degrees of freedom. The second, highly nontrivial step, is to check the positivity in this subspace. This has in se nothing to 
do with having the BRST at one's disposal or not. From parts of the existing literature, sometimes unattentive readers might 
get the impression that having a nilpotent BRST charge is sufficient to ensure unitary of the underlying field theory, as long as 
the BRST is a symmetry. This is actually never true, even at the perturbative level. Indeed, one always needs to check that the 
degrees of freedom singled out by the BRST cohomology do have positive norm. This is, by the way, clearly mentioned in the 
classic work [33], see the discussion below eq. (3.6) on p. 25. An explicit counterexample of a gauge theory with a nilpotent 
BRST invariance which is not unitary in terms of its (asymptotic) elementary degrees of freedom can be found in e.g. [43]. 
Things get even more complicated in the context of a confining gauge theory: what if one cannot define the asymptotic states of 
the elementary excitations when the latter are confined? According to our understanding, one then needs to access the spectral 
properties of the two point correlation functions of composite operators, constructed from their confined constituents. It is now 
clear that the BRST symmetry itself does not provide any answer to such inherently nonperturbative spectral question. In fact, 
it appears that very little is concretely known about the spectral properties of glueballs or hadrons in general. For example, 
lattice QCD offers a powerful way to access gauge theories nonperturbatively, but from a Monte Carlo simulation in Euclidean 
space time at momentum p 2 > 0, it is not straightforward to extract anything on the behaviour of e.g. two-points functions at 
complex momenta. Also, most functional approaches to the bound state problem, viz. a combination of Dyson-Schwinger and 
Bethe-Salpeter equations, are usually not focused on obtaining information in the whole complex momentum plane. Anyhow, 
given the nontrivial analytical structure of the input Green functions of gluons, ghosts and quarks, it is to be expected that it 
will not be a trivial matter to assure the correct spectral properties of their bound states. For example, propagators with complex 
conjugate poles lead to cuts in the complex plane when these poles are suitably combined, see e.g. [12] for an exploratory study 
in a Gribov-like toy model. 

In order to avoid dealing directly with the spectral representation of the correlation functions of composite operators, an 
alternative approach would be to construct, directly from the QCD action, the effective Hamiltonian for e.g. the glueball degrees 
of freedom, and to carry out a canonical quantization with the latter Hamiltonian. Unfortunately, as far as we know, this is far 
beyond the current capabilities. In general, we are talking here about particles with higher spin, notoriously hard to handle in 
a self-consistent relativistic quantum field theory way. The next best thing would be a non-relativistic Hamiltonian description. 
But also here, who knows how to construct even this ab initio from the QCD action? 

To end, we believe a lot of work needs to be carried out before a full satisfactory understanding of the complicated issue of 
unitarity in a confining gauge theory will become available. With this letter, we hope to have provided insight into at least one 
face of this problem: even with a spontaneous breaking of the BRST as is the case in the nonperturbative Gribov-Zwanziger 
Landau gauge quantization, a consistent definition of a time invariant operator subspace can still be given, very much alike as in 
the case of unbroken BRST. As a consequence, the unitarity question in both unbroken and broken case boils down to the same 
difficulty: to establish whether the composite operator Green functions do display the correct analytical properties to represent 
true physical degrees of freedom. 
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